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Abstract. An investigation of the kinematics of a plasma stream rotating in the pulsar magnetosphere is presented. 
On the basis of an exact set of equations describing the behavior of the plasma stream, the increment of the 
instability is obtained, and the possible relevance of this approach for the understanding of the pulsar rotation 
energy pumping mechanism is discussed. 
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1. introduction 

One of the important stages in the development of pul- 
sar radiation theory was the discovery that the rotation 
energy could transform into energy in the electrostatic 
field; this possibility stimulated the modeling of pulsar 
magnetospheres (see Ref. IK ■ I n a typical model, the 
electrostatic field E , generated over the pulsar surface, 
has a finite projection along the dipole magnetic field Bo- 
The resulting electric force, eEg exceeds the gravitational 
force near the neutron star surface, and ends up uprooting 
charged particles from the surface layer. These particles 
are accelerated in the electric field and radiate 7 photons 
because of the curvature radiation in the dipolar field. 
For photon energies greater than twice the electron rest 
mass ( e 7 > 2m e c 2 ), the electron-positron pair production 
(7 + Bo — > e + + e~ + 7') in the ambience of the magnetic 
field becomes possible. The particles, so produced, repeat 
the cycle; they arc accelerated, generate 7 radiation which 
creates more electron-positron pairs. This cascading leads 
to the build-up of a large flux of relativistic e + e~ pairs 
; the process continues till the pair plasma screens the 
electric field Eo as shown in Ref. El 

The e + ,e~ population in the pulsar atmosphere may 
be conveniently divided into three principal components: 
1) the basic plasma mass (Bulk) with concentration n p i, 
and a Lorentz factor j p i, 2) a tail with concentration n*, 
and Lorentz factor 7 t , and 3) the remnant of the pri- 
mary electron or ion beam (most likely electrons). Note 
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that cascading due to pair production is possible only in 
a highly relativistic scenario, because for nonrelativistic 
or even mildly relativistic velocities, the pair annihilation 
dominates pair production (see, for example, Ref. |SJ|SJ). 
The reader is referred to Ref. 013 for a guide to processes 
that define and control the behavior of radiation in the 
pulsar magnetosphere. 

The essence of this "standard model" is that the pul- 
sar magnetosphere is permeated with a multi-component 
plasma. It is generally assumed that the energy con- 
tents of the three major components are of the same or- 
der: ripijpi ~ n t "ft ~ nb'Yb- Theoretical estimates for the 
density and the energy of the primary component are: 
rib ~ (10 11 — 10 12 )cm~ 3 ( also called the Goldreich Julian 
concentration) and 7& w 10 6 — 10 7 . The origin of the pul- 
sar radiation is supposed to be in the magnetosphere, and 
hence the energy contained in the e + e~ plasma must be 
enough to account for the observable radiation. 

In these pulsar models, the measure of the region over 
which the electric field is nonzero is in great significance. 
The distance between the star surface and the region 
where it is screened out is called the vacuum gap, and 
may be estimated to be (10 3 - 10 4 )cm for E Q w 10 7 G(see 
for example Ref. |UJ). Unfortunately the particle energy 
inventory accumulated within the gap is not sufficient to 
explain the observed visible radiation. Several mechanisms 
have been invoked to increase the gap size. The interme- 
diate formation of positronium (electron-positron bound 
state) which would hold back the decay, and lead to an 
increase in the gap was proposed in Ref. 0]] According to 
Ref. Fl, the pulsar magnetic field lines, supposedly curved 
towards rotation, must be screwed even more; the twisted 
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magnetic field lines will be rectified, and as a result will 
enlarge the gap. A different mode, the PFF (pair forma- 
tion front) mechanism was introduced instead of the gap. 
Taking the heating of the stellar surface and the electron 
thermoemission into account leads to nonstationarity, and 
explains nonstationary behavior of pulsars (see Ref. fT2 j ) . 

In a series of papers, Muslimov and Tsigan have at- 
tempted to solve the gap problem via the general rela- 
tivity route. Realizing that in the vicinity of the rotating 
neutron star, the space-time is slightly curved, they work 
out the creation of the electric field in Kerr metric (see 
Ref. 13). As a result, the gap size is somewhat enlarged 
but not enough to make a difference. 

It seems that some additional source of energy or some 
new mechanism will be necessary to surmount the prob- 
lem arising out of the "insufficiency" of the energy con- 
tent in the vacuum gap. For instance, it seems perfectly 
plausible and possible that one could draw upon the pul- 
sar rotation energy to augment the energy content of the 
magnetospheric e + e~ plasma; the rotation energy, for ex- 
ample, could be transformed into energy associated with 
oscillations in the e + e~ plasma far from the pulsar sur- 
face. In this region the complicating effects of unremovable 
gravitation (having to use, the Kerr metric, for example) 
will not exist. This paper is an attempt to formulate and 
investigate the problem associated with the parametric 
pumping of plasma oscillations in the pulsar magneto- 
sphere taking into consideration the complicated nature 
of the multicomponent electron-positron substrate. 

Let us review a standard system to examine the plau- 
sibility of the proposed mechanism. For the Crab Nebula, 
the nebula radiation can be surely sustained by the ro- 
tational energy of the pulsar PSR 0531 located in the 
vicinity of the nebular center. The power of the nebula 
radiation exceeds the pulsar radiation power by two or- 
ders of magnitude (it is approximately 2 • 10 38 erg/sec). 
The only source capable of providing such a prodigious 
power is the slowdown of the pulsar rotation: W « JQf2, 
where / is pulsar's moment of inertia with mass of order 
(1.5 — 2.5)M (Mq is solar mass), f2 and Cl are star's an- 
gular velocity and angular acceleration respectively. This 
power is equal to 5 • 10 38 erg/sec. The rotation energy loss 
rate is described by the ratio W /W ~ 2VL/Q = 2P/P, 
where P — is the neutron star's rotation period, 

and P and P are measurable parameters for pulsars. Ratio 
P/P for different pulsars ranges from 10~ 11 sec _1 (PSR 
0531) to lO-^sec" 1 (PSR 1952+29). Thus, for ordinary 
plasma oscillations to be pumped by the pulsar rotation, 
their growth rate must not exceed P/P, the measurable 
rate of rotation energy loss. 

A hydrodynamic approximation will be used to study 
the problem of rotation induced wave generation in plas- 
mas. For relative simplicity one will assume that the e + e~ 
plasma has two distinct energy ranges: the lower range 
(still highly relativistic) with n p i, and 7 P ;, and the beam 
with rib, and 7&- These flows propagate along the rotating 
monopole like magnetic field lines. 



2. The main consideration 

It is supposed that for distances less than the radius of 
curvature of the field line, the magnetic field is monopole 
like (in this approximation the magnetic field lines may 
be supposed to be rectilinear). The problem of the mo- 
tion of charged particles in the pulsar magnetosphere can 
be considered in the local inertial frame of the observers, 
who measure the physical quantities in their immediate 
vicinity. They are called the Zero Angular Momentum 
Observers (ZAMOs). Naturally, the proper time of the 
observer riding the particle, is different from the proper 
time of ZAMOs. 

Transformation from the inertial frame to the frame 
connected with the pulsar magnetosphere can be done in 
the following way: 

t = t',<p = ip',r = r',z = 0, (1) 

then the interval in the corotating frame will have the 
form: 

ds 2 = - (1 - fiV) dt 2 - dr 2 (2) 

where VL is angular velocity of rotation. 

Note that the lapse function a = \fl — fl 2 r 2 (c = 1) 
not only connects the proper time of ZAMOs with the 
universal time dr = adt, but also gives a gravitational 
potential: 

Va 

g = ■ (3) 

a 

The equation of motion of the particle in the ZAMOs 
frame is expressed as: 

^= 7 g+-(E+[VB]) (4) 
dr m 

where 7 = (1 — V 2 )~ 1//2 is the Lorentz-factor and V = 
dr/dr is the velocity of the particle determined in the so- 
called 1+1 formalism, and p — » p/m is the dimensionless 
momentum. 

As is shown in Ref. 1 141 the transition from the particle 
equation of motion to the Euler equation for fluid dynam- 
ics in the 1 + 1 formalism, may be easily fulfilled if one 
changes d/dr in Eq. (4) by l/(adt) + (VV). The result- 
ing equation describing the stream motion (neglecting the 
stream pressure) takes the following form: 

+ (VV)p = -7— + - (E + [VB]) (5) 

a at am 

where V and p are now hydrodynamic velocity and mo- 
mentum respectively. In order to rewrite this equation in 
the inertial frame, let us note that the ZAMOs momen- 
tum coincides with the momentum in the inertial frame. 
In fact, from the definitions p = 7V, 7 = arf ', and 
V = dr/dr (prime refers to quantities in the inertial 
frame), one can easily find that p = p'. In the inertial 
frame, then Eq. (5) converts to (omitting primes for all 
quantities) : 
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dpi 

~dt 



(viV)pi - - 7 aVa + — (E + [v ; B]) , (6) 

m 

i = b,e,p 



where 6, e and p denote the beam, electron and positron 
components respectively. In Eq.(6), the force F = -7«Va 
is the analog of the centrifugal force. 

Adding the continuity and the Poisson equations: 



^i+V(n 4Vi ) = 0, 



(7) 

(8) 



VE = 47re(n e — n p + rif,) 

As we are interested in the evolution of fluctuations, 
one can look for solutions of Eqs. (6)-(8) in the framework 
of a perturbation theory- an expansion in which terms 
like E l /mwy (the small parameter in the approximation 
of weak turbulence for the plasma) are small: 



E = E 



E 



(9a) 



B = B° + B 1 



(96) 

Pi = P?+Pf + - (9c) 

where E°, B° and p? are the leading terms, and E 1 , B 1 
and p?" constitute the perturbations. In the zeroth ap- 
proximation (taking into consideration the fact that the 
ejected particles not only move along the radius, but also 
corotate with the pulsar magnetosphere because of the 
frozen-in condition , Eo + [voBo] = 0) Eq. (6) will be 
reduced to the form (see for example Ref. 1151 IT?))) : 



d r 



i - n 2 




(10) 



where we have neglected the term (voV)po (taking this 
term into account is a separate problem, and is not exam- 
ined within the framework of this paper). As is shown in 
Ref. El Eq. (10) allows an exact solution for particular 
initial conditions r(to = 0) = 0, V(to = 0) = Vq: 



r(t) = 



V Sn{Qt I m) 



(11) 



dn(£lt | m) 

where Sn and dn are Jacobian elliptical functions, the sine 
and the modulus respectively and fh = 1 — Vq. 

Using the following properties of the mentioned 
Jacobian elliptical functions Sn(x | 0) = sin(i) and 
dn(x | 0) = 1 (see Ref. IT7jl . one may easily reduce Eq. (11) 
for the ultra-relativistic regime (Vq — ► 1) relevant to this 
paper: 



r(t) = sin fit. 



(12) 



With this known asymptotic solution, the first order 
Eq. (6) reads: 



M + (Tr v)p? 



-E 1 



(13a) 



Fj = n 2 rv°pl, 

v i = Vio cos fit 



(136) 
(13c) 



where, in addition to Eq. (12), we have used the Lorentz- 
factor expansion in the small parameter pj/p®: 7 



1 + pf (1 + P?pf/(1 + P? )j 
In a similar fashion, the linearized Eq. (7) becomes: 



^+div(n°v 1 1 ) + divK 1 v°)-0. (14) 

The electron-positron continuity equations may be 
combined to obtain the evolution equation for the effective 
charge density n p i = n e — n p : 



dn 



P i 



Of 



divKv^) + div(n^v^) = 



(15) 



where Vpj = — vj and it has been assumed that v° = 
v° p = v°, and n° e = n° p . 

We seek here a solution in which the density pertur- 
bations have no spatial dependence so that the last terms 
of Eqs. (14) and (15) are identically zero, let us choose n\ 
(now i = pl,b) to have the form 



1) = Ar-e~^ sinnt . 



(16) 



Then, if the rest of the perturbed quantities are al- 
lowed the spatial dependence p*(r;t) = p*(k; t)e lkr , one 
obtains from Eqs. (14)-(15): 



ni dt 



(17) 



where 70,: is the initial Lorentz factor, and Ri = ^ sin £lt, 
and the relation vj = p]/^ {p pl — pl~P p ) that is satisfied 
for fit ~ Sl/ui <C 1 has been used. 

Let us now go back to Eqs. (13) and write them sepa- 
rately for the two components- the plasma and the beam 



dpi 

dt 



+ Kv) p i 



F£ + -E J 

m 



(18a) 



+ (v»P^ 1 =F^ 1 + 2|-E 1 . (186) 

Repeating the procedure applied to Eq. (18) and tak- 
ing Eqs. (12) and (17) into consideration, one can find 
that the beam and the plasma components evolve as: 



d 2 N b 
dt 2 



l b „ikR b 



kEt, 



d 2 N, 



pi 



dt 2 



= -2i 



m 1p0 



e ikR ^kE x . 



(19a) 



(196) 



Combining Eqs. (19a) and (19b) eliminates the electric 
field to yield 



d 2 N h 



dt 2 



(20) 
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Equation (20) is a rather complicated non-autonomous 
equation in time. To solve it, one can take its Fourier time 
transform (restoring the speed of light) along with that 
of the Poisson Eq.(8) [cf. Appendixes A,B Eqs. (A. 2) and 
(B.2)] with the electric field eliminated to arrive at the 
coupled system 



where sq = [Lu p i/Q,j^ p f]. In addition to the real root, the 
dispersion relation of Eq. (26) allows the complex conju- 
gate pair 



Ai, 



1 iV% 

2 2 



(27a) 



uu 2 N b (uu) 



3 

^ELJ2^ + sn ) 2j s(a)N pl (u + sn), (21) 

Zn pHob s 



-^L\N p (cj) = -fJ2 J s(a)Nb(u-sn) (22) 
7o P i / 7o P / s 



M 



3 3/2 s 



Jl(a) 



1/3 



(27b) 



with comparable real and imaginary parts. The root with 
the positive imaginary part implies the instability that we 
were seeking. 



where a = kc/2QjQ pl , and uj p i — ^J8nn pl e 2 /m. Naturally 

one can see the appearance of convolution sums on the 
right hand sides of both equations. 

Substituting iV& from Eqs. (21) into (22), one can find: 



7o P / 



^£j s (a)J z (a) 



7o& 



lu — sQ 



N pl (u - (s - 

(23) 

where Lu b — y^Sim^e 2 /m. One could try to find the na- 
ture of the time evolution from Eq.(23), but it seems 
to be a little better to go back to n p i given by N p i = 
n pl e lkRpl . Carrying out the algebra given in the ap- 
pendix[cf. Appendix C, Eq. (C.2)], one may derive 



N pl = 



2J J s (b)n p i(uj - sti) 



(24) 



where b = kc/fl. Substituting Eq.(24 ) into Eq. (23), one 
finally obtains the rather complicated dispersion: 




Fig. 1. Graph of log(S(R/ R cy i)). Set of parameters is fol- 
lowing: n ob fa 10 7 cm~ 3 , 7oh « 10 6 , ~f 0pl ps 10, n ps 190 Hz, 
r pa 10 6 cm. 



LU 



J s (b)n p i(oj - sQ) 



^ Ji(a)J J+ i +m (a)J m (b) 



r6i 



jlm 



lu + (j + m)f2 



n p i (oj+jtt). 
(25) 

To extract some sense out of the above result, one 
can explore the the dispersion relation near the resonant 
condition, lu 2 ps ou 2 i/7o P i- Near the resonance, the basic 
contribution to the sum ^ s J s (b)n p i(uj — sfi) comes from 
lu ps sq^I. Similarly the right hand side of (25) is reduced 
to a single term corresponding to lu pa — Jo^, oj ~ 
and m = m = -j = l Q = s . Rewriting uj 2 - oj 2 J^l pl 



3/2 



as 2u p iA/i Qpl 



3/2 



(where A = oj - oj p i/i 0pl 
Eq. (25) the simple cubic equation: 



u 2 uj p i 



/2 j 

3/2 so 



a) 



), then, reduces 



(26) 



3. estimates 

The instability turns out to be rather strong. For typi- 
cal pulsar parameters (that of the Crab nebula) nob ~ 
10 7 cm" 3 , job « 10 6 , 7o P / w 10, fl w 190Hz, consider- 
ing the longitudinal waves, one can plot the graph of the 
growth rate 5 (S = \fiMj2 see Eq. (27)) as a function 
of the radial distance (normalized by the light cylinder 
radius). For simplicity let us examine longitudinal waves 
in a coupling point. In this case all three modes of a cold 
plasma (o, x and Alfven modes) are indistinguishable, then 
one can estimate the wave number according to the fol- 
lowing approximate formula kc ~ ^pi/lopf ( see R- ei -EJ)- 
Using mentioned pulsar parameters, from Eq. (27a) the 
growth rate is estimated. In Fig.l one may see that 6 
is sensitive to the radial distance, in the range 0.75 < 



R/Rcyi < 1 it increases from 



io- 



to 



io- 



(here 



i? cy i is the light cylinder radius, thus the radius, where the 
rotation velocity equals to the speed of light). It is clear 
that, S becomes unreasonably large for R/R cy i > 0.85, 



Osmanov et al.: Parametric mechanism of the rotation energy pumping by a relativistic plasma. 



5 



which means that, the linear assumption will be grossly 
violated long before these distances are reached. The insta- 
bility will have a linear growth over certain distances and 
then nonlinear saturation will set in after which other non- 
linear mechanisms may take place in the energy pumping 
phenomenon. In estimating the growth one has invoked 
no p ijopl ~ «ofc7ob, and used the fact that the density goes 
as n w n (r/R) 3 (where r w 10 6 cm is a radius of the 
neutron star) because the magnetic filed is monopole like. 

4. conclusions 

The purpose of this paper was to explore the possibility 
of pumping the rotational energy of the pulsar into the 
plasma. Considering a highly idealized system, the linear 
instability caused by rotation in a two component rela- 
tivistic plasma embedded in a uniform magnetic field was 
examined. By using the hydrodynamic equations of mo- 
tion, the continuity and the Poisson equation, it has been 
shown that the plasma waves can grow on the rotation 
energy with rather high growth rates. In fact the pertur- 
bation growth rate for distances R/R cy i > 0.85 is quite 
high in comparison with the supposed instability rates we 
were seeking. If this mechanism is, indeed, operational, 
then the only consistent scenario is that the linear stage 
with this growth rate is very short, and nonlinearities are 
turned in soon enough to considerably reduce the growth. 
Thus the need for a nonlinear theory is immediately and 
strongly indicated. Sooner or later we should consider this 
particular nonlinear effect, which will comprise one more 
step closer to the real scenario. There are two other short- 
comings of this effort: a) the straight magnetic field lines 
has been examined, whereas real profiles are curved and 
b) only electrostatic waves have been considered. The pre- 
liminary results, however, unambiguously show that the 
energy content of the magneto spheric plasma can grow 
at the expense of the stellar rotational energy. 
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Appendix A: Derivation of Eq. (21) 

Note that: 



where 
then 



N f (t) = J du'e-^Nfiu') 
f = b,pl 

dtdcj'e lu ' t (-iuj') 2 e- luJ ' t N f (cj') 



(A.la) 



= - J J dtdLu'e lt ^'\uj') 2 N f {uj') = 
= -2tt J du)'(J) 2 5(w-u')N f (uj') = 



= -2ituj 2 N } (w) 



(A.lb) 



where the following representation of the delta function 
has been examined: 



5(x) = -!- / dke lkx 



tipi — Kb = — sin i It : 



(Ale) 



fi 



C 1 1 



sin fit 



■ — 2— sin fit, (A.ld) 



then 



g—ik(R p i—Rb) ^ ^iasinQt 

kc 



a 



(Ale) 
(Al/) 



We have considered the approximate expression: Voi ~ 
1 — l/27o i; and the following observable fact: -f b 3> 7o P ;- 
As one can see, the speed of light has been restored again. 

By using the following identity: 



±ix sin Clt 



,±isOt 



(Al<?) 



one may easily transform e lk ( R p R <>) into its Fourier 
mode: 



if- 



= -Y^J s (a) J J dtdLu'(Lu') 2 Np(Lu'y t ^ +sn -"') = 
= -2nJ2 J s(a) I duj'(J) 2 Np(uj')5(u + sn-J) = 

s J 

= -2irJ2 J s (a)(uj + stl) 2 N p (lu + sty. (A.lh) 



Combining Eqs. (A. lb) and (A.lh), taking into con- 
sideration Eq. (20), one finally will obtain: 



3 

u; 2 N b (u;) = pfi^ ^(c + sO) 2 J s (a)N pl (cj + S n). (A.2) 

ln pHob s 
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Appendix B: Derivation of Eq. (22) = _ 2?r ^ Jg{Af) J duJ ' nf{uj r )5{uj + sQ 

One can easily transform Eq. (8) into the following form: s 

kE 1 =4*e(N b {t)e- ikR >+N pl (t)e- ikR "). (B.la) = -2^J s {A f )n){u - s^l) 

Taking Eq. (19b) into account, one can have: Hence we will have: 

= (N A t) + N b ( t) e^-«»A . N 'M = S WfWu, ~ -A) 

dt 2 m % p i 

(B.lb) 

For Fourier expansion of the left-hand side of 
Eq. (B.lb) one can analogously go to Eq. (A. lb) write 
down: 

/ J ' du'e-^'Nniu') = -2ituj 2 N p1 (u) 

(B-lc) 

while, for the second term in a bracket, one finds: 
^J s (a) J J dtdLu'e u ^- sn -"'^N b (Lu') = 

= -2nJ2M a ) ( duj' N b {uj')5{uj - stt - J) = 

s J 

= -2irJ2 J s (a)N b (uj - sty. {Bid) 

s 

Combining Eqs. (B.lc), and (B.ld) we will obtain 
Eq. (22): 

u? - -f) N pl (uj) = ^^J s (a)Ar fe (u;-sQ) (B.2a) 
"topi J lopi s 

where 
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Appendix C: Derivation of Eq. (24) 

From Eq. (16) we have: 

N f (t)=n)(t)e ikR f. (C.la) 

Taking into consideration Eq. (A.lg), and a condition 
7o ~> 1 one finds: 

N f (t)=J2j s (Af)e isQt n}(t) (C.lb) 

S 

where 

A t = £ (C.lc) 

Fourier analysis of the right-hand side of the Eq. (C.lb) 
gives following: 

J^MAf) J J dtduj'e it ^ +sn - u '' ) nf{uj') = 
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